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Abstract 

It is well known that any amount of energy injected in a harmonic oscillator which 
is resonant and weakly coupled with a second harmonic oscillator, tunnels back and 
forth between these two oscillators. When the two oscillators are anharmonic, the 
amplitude dependence of their frequencies breaks, in general, any eventual initial 
resonance so that no substantial energy transfer occurs unless, exceptionally, an al- 
most perfect resonance persists. This paper considers this interesting situation more 
generally between two discrete breathers belonging to two weakly coupled nonlinear 
systems, finite or infinite. A specific amount of energy injected as a discrete breather 
in a nonlinear system (donor) which is weakly coupled to another nonlinear system 
(acceptor) sustaining another discrete breather, might be totally transferred and 
oscillate back and forth between these donor and acceptor breathers. The condi- 
tion is that a certain well defined detuning function is bounded from above and 
below by two coupling functions. This targeted energy transfer is selective, i.e., it 
only occurs for an initial energy close to a specific value. The explicit calculation 
of these functions in complex models with numerical techniques developed earlier 
for discrete breathers, allows one to detect the existence of possible targeted en- 
ergy transfer, between which breathers, and at which energy. It should also help 
for designing models having desired targeted energy transfer properties at will. We 
also show how extra linear resonances could make the energy transfer incomplete 
and irreversible. Future developments of the theory will be able to describe more 
spectacular effects, such as targeted energy transfer cascades and avalanches, and 
energy funnels. Besides rather short term applications for artificially built devices, 
this theory might provide an essential clue for understanding puzzling problems of 
energy kinetics in real materials, chemistry, and bioenergetics. 
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1 Introduction: a brief review of discrete breather theory 

It is a great pleasure to contribute to the volume dedicated to Professor E. N. 
Economou with work that attempts to bridge the gap between localization and 
propagation in complex discrete systems. In the long-standing contributions 
of Prof. Economou, localization is induced via quenched disorder while prop- 
agation is a result of either higher dimensions or nonlinearity. In our present 
work, coherent energy transfer occurs as a result of judicious interlocking of 
disorder with nonlinearity and an appropriate exploitation of the concept of 
nonlinear resonance. Resonance is the basic principle for energy propagation 
in spatially periodic linear systems. The simplest example in textbooks is ob- 
tained for two coupled identical harmonic oscillators. In this case, it is well 
known that any amount of energy deposited on one of the two oscillators tun- 
nels back and forth between these two oscillators with a frequency proportional 
to the coupling. On the other side, in order to keep the energy localized in a 
given harmonic oscillator, the resonance between this oscillator and the other 
harmonic oscillators must be broken. This situation can be reached in a har- 
monic infinite system for an impurity mode or, more generally, by breaking 
translational invariance, for instance through strong enough disorder or by 
quasiperiodic modulations. Roughly speaking, finding a harmonic oscillator 
that is well resonant within the frequency interval Sou with another oscilla- 
tor generally requires to go far from the initial one, typically at a distance 
5uj~ 1 I d in a model at dimension D. But since the effective coupling between 
these two almost resonant oscillators drops exponentially as a function of their 
distance, the resonance is generally not sufficient for allowing a substantial en- 
ergy transfer. This effect blocks energy propagation and causes the Anderson 
localization of the eigenmodes [1,2]. 

When the system becomes nonlinear, this picture is drastically changed be- 
cause the frequency of an anharmonic oscillator depends on its amplitude. 
Concerning the time-periodic solutions, we have shown that nonlinearity plays 
a "double game" [3-5]. First, by adjusting appropriately the amplitude of 
the nonlinear oscillators, nonlinearity could restore resonances between an- 
harmonic oscillators which were lost in the linearized system. Indeed, small 
nonlinearity is sufficient to generate exact extended solutions in random sys- 
tems while the whole linear spectrum is localized. Although these solutions 
appear spatially intermittent, they can propagate some nonvanishing amount 
of energy (much smaller though than standard uniform plane waves). Their 
existence demonstrates that random systems with a small nonlinearity can 
propagate energy while the same systems without nonlinearity cannot. Sec- 
ond, with different amplitudes for the nonlinear oscillators, nonlinearity can 
also maintain (or magnify) the lack of resonance of a localized vibration with 
the other modes. Thus, the same system can also sustain exact localized non- 
linear modes. Moreover, the existence of these modes may not require random- 
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ness. Mode localization can also be achieved in a nonlinear spatially periodic 
system provided the frequency and its harmonics are nonresonant with the 
extended linear modes (i.e., they do not belong to the linear phonon band). 
These modes are called discrete breathers [6,7]. 

Discrete Breathers (DBs), also called Intrinsic Localized Modes, are spatially 
localized time-periodic solutions of discrete classical nonlinear dynamical sys- 
tems which could be finite or infinite in arbitrary dimension, spatially periodic 
or random, etc. Many early works predated their discovery. For example, po- 
larons are nonlinear objects belonging to the family of DBs introduced long 
ago by Landau [8]. The concept of local modes was independently introduced 
in chemistry also long ago [9-11] for small molecules. In that case, discrete- 
ness was implicitly necessary. Incidentally, these solutions are recognized to 
be important for understanding the quantization of molecular vibrations [12]. 

Renewed interest on these nonlinear concepts but within the framework of inte- 
grable models which have exceptional (but sometimes physically misleading) 
properties started during the 1970's. The self-trapping equation introduced 
later in the mid-1980's already revealed some features characteristic of DBs 
but only in special cases (see a recent interesting review by A. Scott [13]). 
Many interesting physical ideas of applications in physics were proposed, but 
were limited by flaws of the underlying theories. They could be reconsidered 
and amended using the more universal concept of DB. 

It was only in 1988 that Sievers and Takeno [6] claimed the existence of DBs 
as long lifetime solutions in general models, finite or infinite, on the basis 
of approximate analytical and numerical calculations (see review [14]). Actu- 
ally, their claims were definitely confirmed later by a rigorous proof of their 
existence as exact solutions in Klein-Gordon systems under rather general con- 
ditions of non resonance [15]. Moreover, DBs were proven to be linearly stable 
and robust to model perturbations (i.e., not restricted to specific models). 
Subsequently, the scope of validity of this proof has been extended to more 
general models with optical or acoustic phonons and with arbitrary complex- 
ity in any dimension [16-20]. This general theory also brought new methods 
for systematic and accurate numerical investigations [21]. 

The existence of DBs essentially requires discreteness and nonlinearity (ran- 
domness is not necessary despite the fact that DBs could persist in its pres- 
ence). Thus, these conditions fully take into account two essential character- 
istics of real matter. It is discrete, because it is made of atoms, and nonlinear, 
because the interatomic potentials are nonlinear. DBs may exist generically 
as exact solutions when the frequency of a local mode and all of its harmonic 
frequencies either belong to phonon gaps or are above the whole phonon spec- 
trum. Otherwise, these local modes would radiate away their energy by the 
fundamental frequency or some harmonics and thus survive only over a finite 
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time. However, if the radiation occurs by high order harmonics, the lifetime 
of this local mode at K may be long but nevertheless finite. Thus, it is clear 
that DB existence as exact solution crucially depends on the gap structure of 
the phonon spectrum. An exact DB solution at a given frequency (or action) 
may exist in one model but may disappear in another very similar model just 
because some high order harmonic becomes resonant with phonons. Another 
essential property is that DBs come as one parameter family which can be 
parametrized by their frequency (or, better, their action which can be used 
for semiclassical quantization) Q 

DBs may have a lot of other amazing properties not fully explored yet but 
which depend on the considered models. They scatter elastically linear phonons 
(without energy exchange) like static impurities. They may filter phonon noise 
coloring the spectrum [22]. Inelastic interactions with small amplitude phonons 
(generally) occur at higher order. Then, DBs may slowly decay while in other 
situations unexpectedly they may grow [23] by pumping energy from the back- 
ground. They are generally pinned to the lattice but may become highly mobile 
in some special situations [24] and then they contribute to energy propagation 
instead of trapping it. They can also react by fission or fusion [25]. 

Because they are exact solutions at K, and the interaction between DBs 
drops very fast as a function of their distance, at low temperature the only 
physical way for a quasi-isolated DB to decay is by the small interaction with 
its own phonon background. Actually, the real spectrum of this phonon back- 
ground is essential for determining the true DB lifetime. Phenomenological 
interactions modeled for example by a white Langevin noise and damping (or, 
equivalently, a coupling with a large collection of harmonic oscillators with 
a uniform structureless spectrum), systematically destroy the DB as an ex- 
act solution simply because the phonon spectrum extends to infinity without 
gaps. Thus, the lifetime obtained by such approaches widely used in physics, 
is nothing but the result of the assumptions that are made and might be 
physically misleading. 

Consequently, when the temperature is sufficiently low (compared to the DB 
energy) and the DBs sufficiently far apart, energy trapping by DBs may 
be substantial and persist over extremely long time in apparent violation 
of the standard Fourier law. Therefore, the most important feature charac- 
terizing (pinned) DBs (and distinguishing them from ordinary phonons and 
local modes), is that in appropriate conditions where they can be created (by 

1 The DBs we discuss here are extraband DBs, to be distinguished from intraband 
DBs defined in Refs.[3-5]. Intraband DBs may only exist when the phonon modes 
are spatially localized in quasicontinuation of the Anderson modes. The essential 
difference is that the frequencies (or harmonics) of intraband DBs lie inside the 
phonon spectrum but then these frequencies are restricted to fat Cantor sets. As a 
consequence, the intraband DBs do not form continuous families versus frequency. 
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external excitations or by other processes breaking initially the thermal equi- 
librium), they may spontaneously appear and persist as spatially localized 
dynamical structures with unusually long lifetime. 

Numerical simulations clearly demonstrate that in systems where linearly sta- 
ble DBs exist, they show up spontaneously under thermal shocks [26,27] or 
when starting from initial states far from thermodynamic equilibrium [28]. 
They trap energy over very long time and generate unusual energy relax- 
ation (stretched exponentials). This long lifetime of DBs should be relevant 
for understanding many anomalous energy relaxations especially in disordered 
matter (glasses, polymers, biopolymers) which up to day do not have yet truly 
satisfactory interpretations. 

Up to now, experimentalists did not pay much attention to this new paradigm 
in physics. One of the reasons might have been that DBs could be easily con- 
fused with early concepts inherited from the theory of integrable systems such 
as (nondiscrete) breathers and solitons Q(perhaps because of the terminology). 
Otherwise, the localization or the anomalous relaxation of energy observed in 
real experiments, are usually interpreted by randomness, metastable states 
etc. 

In spite of the general belief, some experimentalists recently observed DBs with 
different techniques and in different systems, for example in coupled nonlinear 
optical wave guides [29,30], in coupled arrays of Josephson junctions [31,32], 
in some nonlinear materials [33], in magnetic systems [34] and even possibly 
in myoglobin [35]. 

After this brief review we now come to the main aim of this paper which is 
to present a novel property termed targeted DB energy transfer that enables 
possible perfect transfer of energy between two DBs. This idea was already 
briefly suggested in [36] and considered in term of phase dynamics in [37]. We 
give here a precise analytic treatment of this effect through a generalized dimer 
model. We should point out that our results do not agree with the scenario 
imagined in [37]. Although, this phenomenon can only exist in appropriate 
models and for appropriate DBs, it could become the key for understanding 
puzzling phenomena of well-focused energy transfer in physics and especially 
in biomolecules. 



2 They only exist in very special integrable models and, moreover, their properties 
are fragile to perturbations in contrast with the highly robust DBs. 
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2 Generalized Dimer Model for the Energy Transfer between two 
Discrete Breathers 

The energy tunneling between two resonant harmonic oscillators is described 
through the linear superposition of the symmetric and antisymmetric modes. 
This is of course not possible when the oscillators are (strongly) anharmonic. 
We show in this section that we can model the energy transfer in this case 
by an integrable generalized dimer model. We shall study more generally the 
energy transfer between two DBs so that the energy transfer between two 
oscillators becomes a particular case. 

We consider an arbitrary anharmonic system (D) (finite or infinite) called 
donor molecule described by some set of variables and conjugate variables 
{pf i iP} an d a Hamiltonian Hoiipf, <?f})- We assume that this molecule can 
sustain DBs and select one of them called donor DB. 

As we know, DBs come as continuous family parametrized for example through 
their action (which is Id = J J2iPidqi/27i integrated on the loop representing 
the DB in the phase space). Then, DB's energy Hd{Id) becomes a function 
of action and its frequency is u b = dH D /dl£). We also assume that this DB 
family is gapless, i.e., Hd{Id) — > for Id — > when the DB is vanishing at 
zero energy^. 

We consider a second arbitrary anharmonic system (A) (finite or infinite) 
called acceptor molecule with properties similar to the donor molecule. We 
select on the acceptor molecule a family of acceptor DB with action I a and 
energy Ha(Ia)- As with the donor system, we assume that acceptor DB is also 
gapless. 

Subsequently, we introduce a weak coupling between the donor molecule and 
the acceptor molecule. These two systems could physically represent two in- 
dependent molecules or two linked parts of the same macromolecule relatively 
far from one another. The weak coupling interaction could either come from 
contact (hydrogen bonds, van der Waals forces), from (screened) Coulomb 
interactions, could be mediated by the solvent, or come directly from the in- 
termediate link of the macromolecule if any, etc. 

We would like to find the conditions for complete transfer of the energy of 
a DB from donor molecule with action It, after some time, to a DB on the 
acceptor molecule. Since the coupling is weak, this transfer will be slow. Then, 
at any intermediate time t, the dynamics of the system will be well described 

3 Note that DBs in space-periodic systems in more than one dimension cannot be 
gapless [38] but we essentially consider here systems where the space periodicity is 
broken by the presence of the donor or acceptor sites. 
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over a time scale long compared to the period by two DBs with action Id on 
the donor and I a on the acceptor. This solution is well represented in terms 
of loop dynamics [25] as an almost invariant closed loop in the phase space 
which projects like DBs on the isolated molecules in the complementary donor 
and the acceptor subspace. Thus, the action of this loop is nothing but the 
sum of the actions of these two isolated DBs. Because of the weak coupling 
between donor and acceptor systems, this loop evolves slowly in the phase 
space. Liouville theorem states that the action of a loop evolving in the phase 
space is time- const ant. Thus, we have 



Another necessary condition is the conservation of the total energy Ej- — 
H D (I D ) + H A (I A ) yielding 



which is a function independent of Id for < Id < It- Differentiating (2) 
with respect to Id readily yields that the frequencies u>d = dHD/dlD of the 
donor DB and uja = dHA/dlA remain identical during the energy transfer, 
i.e., the two DBs remain resonant. 

Actually, this condition (2) does not need to be perfectly fulfilled (since the 
coupling energy should also be involved in the energy conservation) but we 
shall assume that it is almost fulfilled. We wish to describe the coupled donor- 
acceptor system through an effective Hamiltonian expressing the energy trans- 
fer only as a function of the variables Id and I a and their conjugate variables 
Od and 6a- 

However, the restriction to two pairs of variables will be valid only if condition 
(2) is not fulfilled precisely or even approximately by any other DB of the 
acceptor molecule for the same total action It and the same total energy Ex- 
In other words, we assume that the acceptor DB is unique. We shall assume 
also, but only for simplicity, that the donor DB is unique 0. This condition 
is imposed for avoiding more complex situations with bifurcating transfer of 
energy to two or more acceptor DB. It is in principle fulfilled for molecules 
without any special symmetry. 

A priori, the small coupling energy, C(Id, I a-, Qd, Q a) depends on Id = Io + I 
and I a = Io — I, on the angle difference Q = Qd — Q a (conjugate to T) supposed 
to vary slowly because of the above assumptions (weak coupling and almost 
resonance), and on the total angle = 9 D + A (conjugate to I ) which on 

4 Otherwise we should have to involve other donor DBs in the return transfer. This 
problem, which is also interesting and related to energy funneling, will be treated 
elsewhere. 



It = Id + Ia- 



(1) 



Hd(Id) + Ha(It — Id) = E T 



(2) 
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the opposite varies fast at the scale of the DB periods. It is then justified to 
average the coupling over the fast variable 9q (which varies practically linearly 
over relatively long time) and drop the 9q dependence. Next, it is convenient to 
split the coupling energy as a sum C(I , I, 9) = C (I , I) + C[(I , I, 9), where 
Cq{Iq,I) is the average of C(I ,I,9) over 9, while the average of C[(I ,I,9) 
over 9 is zero. 

Then, the effective Hamiltonian regarding the two DBs of the coupled molecules 
can be written in the form 

H(I , I, 9) = H (I , I) + V(I , I, 9) (3) 

where H (I , 1) = H D (I + 1) + H A (I - I) + C (I + 1 , 1 - 1) and V(I , 1,6) = 
C[(Iq + I,Iq — 1,9) which has zero average with respect to 9. Moreover, for 
I = ±Jo, either the acceptor DB or the donor DB vanishes which removes any 
9 dependence. Then, V(I ,±I ,9) = 0. 

Since this Hamiltonian does not depend on 9 , I = —&H/d9 = 0, which 
implies that the total action + I a = 2I is conserved. Because of the aver- 
aging over the fast variable 9 , 9 essentially represents the difference between 
the DB phases on the donor and acceptor molecules. Effective Hamiltonian 
(3) with J a constant parameter is nothing but the Hamiltonian for the DB 
phase (see ref. [37] for a formal theory of phase dynamics). 

For donor and acceptor molecules with explicitly given Hamiltonians and cou- 
pling, it is possible to calculate numerically this effective Hamiltonian and thus 
to predict the possibility of complete targeted transfer. This Hamiltonian can 
be formally written as a generalized Discrete Nonlinear Schrodinger (DNLS) 
dimer model when defining the complex variables ipo = + I e -*( 6l o+f)/2 
at the donor site D and ipA = Vh ~ I e -*( e o~9)/2 ^ ^ e acceptor site A. 
Then, the conservation of the total action I is equivalent to the conservation 
of the total norm \iPd\ 2 + IV'aI 2 - This approach could be easily modified when 
the frequency of one of the two DBs is resonant with a harmonic of the other 
one. We could then analyze more accurately for DBs the problem of Fermi res- 
onance usually considered for quasilinear modes. This approach can be also 
extended when there are more than two weakly coupled nonlinear oscillators. 
We then obtain a generalized DNLS model on a lattice. It is also extendable 
for quantum anharmonic oscillators. We then obtain a quantum generalized 
DNLS model which is nothing but a boson model on a lattice. In this case, 
the assumption of conservation of the classical action becomes conservation 
of the number of bosons and nonlinearity becomes many-body interactions. 
On this basis we can generalize early calculations based on rotating wave ap- 
proximations, which were initially restricted to specific models (e.g. theory of 
the Davydov soliton, see [13]). Expanding the obtained effective Hamiltonians 
at the lowest significant order simply yields standard quartic DNLS models. 
However, we also learn from our approach that these approximations are only 
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valid in the weak coupling limit, i.e., close to an anticontinuous limit, which 
was not clearly stated in the early approaches. It is thus inconsistent to de- 
scribe later these DNLS models from their continuous limit where they are 
integrable. These extended problems will be discussed later in further publi- 
cations. 



3 Targeted energy transfer solutions in the generalized dimer model 

We search for solutions of the dynamical system defined by Hamiltonian (3), 
which correspond to a total energy transfer. These solutions are such that if at 
the initial time all the energy is on the donor molecule, i.e., Jd(0) = 2 Jo ^ 
and Jyi(O) = 0, it will be completely transferred at a later time on the 
acceptor molecules, i.e., Jd(^t) = and Ia^t) = Jl>(0). Then, 1(0) = J and 
J(tr) = —Jo. We say by definition that when this situation occurs, we have 
targeted energy transfer. 

Since for J = J , as well as for I = — J , the angle 9 is undetermined, the 
topology of the phase space defined by the conjugate variables (J, 9) is the 
one of a 2D sphere represented in cylindrical coordinates 0. Any point P of 
this sphere can be represented by the two coordinates which are the angular 
coordinate 9 (longitude) of its projection in the (x, y) plane and J is its z 
coordinate (latitude). The poles of the sphere are the points where J = ±J . 

There are two time invariants for this Hamiltonian which are Jo and the total 
energy E. Then, 

V(I ,I,9) + H (I ,I) = E (4) 

defines contour lines on the sphere (J, 9) which represent the orbits of the 
dynamical systems in the 2D phase space. 

Since these contour lines are closed loops, J(i) is time-periodic. Actually, the 
solution in the original system consisting of the coupled donor and acceptor 
molecules, is a two site DB the period of which varies periodically with the 
same period as I(t). Thus, the global solution is quasiperiodic. 

The dynamical solutions of the system can be classified in two types. Type 
1 corresponds to the orbits on the sphere which are homotopic to zero, i.e., 
that can be continuously shrunk to zero without crossing a pole. For these 
trajectories, 9 oscillates between two bounds. The limit case where 9 does not 

5 The topology of this sphere was also considered in [37] for two weakly coupled 
pendula. 
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oscillate, corresponds to exact DBs. They correspond to extrema of the energy 
(4) on the sphere at fixed J . 



Type 2 corresponds to the orbits which turn around the axis connecting the 
two poles. Then, 9 rotates by 2tt around the sphere at each period of I(t) 0. 

A trajectory corresponding to a targeted energy transfer between the donor 
DB and the acceptor DB is such that along the orbit, I(t) varies between 1$ 
and —Jo- Thus, it is an orbit connecting directly the two poles of the sphere 
{1,0). It is straightforward to see that a necessary and sufficient condition 
for having a targeted transfer solution at action 7 , is that the system only 
accepts type 1 solutions at this total action. 

Let us explicit more carefully, this condition. Since the two poles must be at 
the same energy, we have 



At the uncoupled limit V(I , I, 9) = 0, this condition is fulfilled when (2) 
is fulfilled. If H (I , I) is not constant, the contour lines are just circles at 
constant / and are of type 2 and thus, we do not have targeted energy transfer. 

When the coupling is turned on, targeted energy transfer requires that there 
exists an action It associated with an energy Et such that Eq.(5) is fulfilled 
for I = It/2 and E = Et- Then, for E = Et, this closed orbit on the sphere 
is defined by implicit equation 



Function €t{I) vanishes at the poles for I = ±It/2 and characterizes the 
detuning between the donor and acceptor DBs during the energy transfer. We 
term it detuning energy function. In some sense, — ey(7) extends the concept 
of Peierls-Nabarro energy barrier, which concerns only static excitations, to 
DBs, which are dynamical nontopological excitations. However, note now that 
the barrier in energy might be positive or negative. With no restrictions on 
the choice of the coupling functions V(It/2,I,9), it is rather easy to find 
sophisticated dependence on 9 (but a bit unrealistic) , such that there cannot 
exist any contour line connecting the two poles. Thus, it is convenient to 
require some properties for the 9 dependence which also make the coupling 
physically reasonable especially if this coupling is small. We shall assume that 
it is sine-like and that is it has only one maximum and one minimum per 
period for —It/2 < I < It/2. This coupling function V(It/2,I,9) has two 
(and only two) zeros per period with respect to 9 since its average is zero. More 

6 In the case of a nonresistive and noninductive Josephson junction (JJ), which can 
be described by a dimer model, type I and class II characterizes the two possible 
states of the J J without DC current and with DC current, respectively [39]. 



Hq(Iq, Iq) = H (I , —J ) = E. 



(5) 



V(I T /2, 1, 9) = H (I T /2, ±I T /2) - H (I T /2, 1) = e T (I). 



(6) 
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generally, Eq.(6) determines either two values modulo 2n for 9 or no values. 
A necessary and sufficient condition for having two solutions 9 + (I) and #_(/) 
which are continuous functions of / on the sphere (/, 9) is simply 

V-(I) = minV(I T /2,I,9) < e T (I) < maxV(I T /2, 1, 9) = V + (I) (7) 

9 6 

Functions 0+(I) and 9-(I) determine the two segments of the closed contour 
line which connects the two poles of the sphere (1,9). The corresponding 
trajectory is a trajectory which realizes a targeted energy transfer. We note 
that during the time tr, the phase difference between the donor and acceptor 
DBs 9±(I) does not rotate over the trigonometric circle but varies in some 
interval. 

We have just demonstrated (for a physically reasonable sine-like dependence 
of the coupling on the phase difference) that targeted energy transfer solutions 
do persist when the detuning energy function is not strictly zero provided the 
detuning energy function is bounded by the max and the min of the coupling 
function. When the detuning function is small, this coupling may be also small. 

The targeted energy transfer solution at / = It/2 may be lost when varying 
the model parameters. It suffices that condition (7) is no more fulfilled. At the 
threshold, there is a certain value I b of / where V(I T /2, 1,9) = V±(I). Then, 
9 = 9 b = 9 + (I b ) = 9_(I b ). Consequently, point (I b ,9 b ) is an extremum of the 
energy (4) on the sphere so that I(t) = I b ,9(t) = 9 b is time-independent, which 
in the initial system is time-periodic. They correspond to exact DB solutions 
of the coupled system made of the donor and acceptor molecules. 

Thus, the threshold where targeted transfer disappears, can be precisely asso- 
ciated with the appearance of a new DB solution of the whole system. When 
approaching this threshold, the targeted transfer solution becomes cnoidal, 
i.e., the time of transfer diverges because the energy transfer slows down in 
the vicinity of I b where the new DB solution will appear. 



4 Explicit solutions in the general quartic dimer model 

We now investigate explicitly the targeted energy transfer solutions in exactly 
solvable dimer models which describe quite well the general behavior in the 
weak coupling limit. DNLS Dimer models were studied in many early works 
[40-42] but only in special cases with no special attention to the detuning 
function. Energy transfer solutions were already found but for relatively large 
coupling in cases where the detuning function (6) was not small. It is thus 
instructive to reconsider a more general quartic DNLS dimer model as an 
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illustration of our general theory. We consider a dimer with Hamiltonian 

n = (\xM 4 + ^il^il 2 ) + (\x2\H 4 + ^m 2 ) - A(« 2 * + rM (8) 

where xi> X2, ^1 and U2 are parameters. A is the inter-site coupling chosen 
positive for convenience, (ip^iipi) and (ip^ifa) are pairs of conjugate vari- 
ables. Defining first the new pair of conjugate variables (li,0i) and (^2,^2) 
by ipi = \p^e 1 and ip 2 — y/T^e" 1 ® 2 and next the pair of conjugate variables 
(J = (/i + / 2 )/2, 9 = O1+O2) and (/ = (Ji- J 2 )/2, = 0i-0 2 ), an equivalent 
form of Hamiltonian (8) is 

W( J, 0, J , O ) = #o(/o, /) + V(I , I, 9) (9) 

with 



H (I , I) = i X o(/o + I 2 ) + XV + w„Jo + wl (10) 



\/(/ , /, 6) = -2\yJtf-PcosO (11) 
and the new model parameters 

X0 = Xl+X2, X = Xl~X2, CJ =^l+^2, ^ = CJi-W 2 (12) 

Complete targeted energy transfer is obtained for 

I T = , E T = — \ uj (13) 

X x\x J 

The detuning function (6) and the coupling functions V + (I) and V-(I) in 
Eq.(7) are 

e T (/) = i X o(/|-4/ 2 ) and V + (I) = -V_(I) = \y/l* - AP (14) 



so that inequality (7) is fulfilled for all I when 

Xou 



A > 



4x 



(15) 



Targeted energy transfer is obtained for any nonzero coupling when xo — 0, 
i.e., when the nonlinear coefficients on the two sites are opposite. On the 
contrary, it is never obtained when x — 0, i.e., when the nonlinear coefficients 
on both sites are equal, except when 00 = 0. The case x = 0, = 0, which 
was treated in Ref . [40] , is special because It is undetermined. Then, there are 
exact energy transfer solutions for any I when the coupling is large enough, 
i.e., A > Xoh/A. Exceptionally, energy transfer is not selective in this case. 
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Fig. 1. Cylindrical projections of the contour lines on the sphere (I, 6) of Hamiltonian 
(9) for I = 0.5 (a), I = 1.0 (b) (ideal targeted transfer) , I = 2.0 (c). Model 
parameters are xo = 0, x = 1) <^o = 0, lo = —1, A = 1/2. 

Eq.(4) for Hamiltonian (9) yields the model trajectories which appear as con- 
tour lines on the sphere (J, 6) defined by 



They are shown for several values of the total action I and for a choice of 
the model parameters such that targeted energy transfer occurs for a certain 
initial action I and energy E (cf. Fig.l). An example is also shown when 
targeted energy transfer never occurs at any action I Q and energy E (cf. Fig. 2 




)■ 



Eliminating 9 between Eq.(16) and the Hamilton equation 




(17) 



yields 



\i 2 + v(i) = o 



(18) 
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Fig. 2. Same as fig.l for I = 0.5 (a) and Iq = 2.0 (b) but with parameters xo = 5, 
X = 1, u = , u = -1, A = 1/2. 

where 

V{I) = \ (\ X0 I 2 + (%/o + w)J + (^xo/o 2 + <Vo " ^)) 2 " 2A 2 (/ 2 " ^) 

(19) 



Eq.(18) describes a single nonlinear oscillator at zero energy with unit mass 
in a potential V(I) Q. The solutions of eq.(18) oscillate in an interval deter- 
mined by two consecutive zeros of V(I) and where V(I) is negative. They are 
determined (apart from a time shift) from the knowledge of the two invariants 
E and Jo and the interval of zeros. 

Since in our case V(I) is a fourth degree polynomial, the number of real zeros 
is even. Note also that all the real zeros I a of V(I) (P(I a ) = 0), necessarily 
belong to the interval —Jo < la < h since outside this interval V(I) is the 
sum of two positive terms in Eq.(19) and thus cannot vanish. 

The situation where V(I) has no zeros never occurs, whatever the initial con- 
ditions are, because I 2 has to be positive at time 0. When V(I) has one pair 
of real zeros, there is a unique solution (apart from a time shift) for I(t) (and 
0(t)) which is time-periodic. It corresponds to a quasiperiodic solution of the 
initial DNLS equation since 

O = j^- = ~2A ° cosd + xolo + XI + w . (20) 

Mo ^//n 2 - P 



7 Note that we have to discard the fake time-constant solutions I(t) = I a , where 
I a is some zero of V(I), which are not solutions of the real dynamical system. The 
reason is that the time conservation of H is trivial for time-constant solutions but 
it is not sufficient to insure they fulfill the dynamical equations. 
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When V(I) has two pairs of real zeros, there are two time-periodic solutions 
I(t) that, again, correspond to quasiperiodic solutions of the initial DNLS 
equation. 

Time-periodic solutions (corresponding to DBs) are obtained when I(t) be- 
comes time-constant (the interval of oscillation shrinks to zero), i.e., when 
V(I) gets a degenerate pair of zeros. 

For solutions corresponding to a targeted transfer of energy between site 1 and 
2, I(t) should oscillate between +1$ and —1$ which requires that ±io are zeros 
of V(I). This condition yields Iq = It/2 and E = Et defined by Eqs.(13). 
Otherwise, V(I) must not have any other zeros, yielding for this parameter 
value 




Vil) -^[l 2 - -) (21) 

which implies 



k 



4A X 



< 1 (22) 



equivalent to condition (15). Then, the solution of Eq.(18) for V(I) given by 
eq.(21) is a Jacobi elliptic cosine with modulus < k < 1 



I(t) = -— cn(Av / 2t; k) (23) 
X 

which describes explicitly the targeted energy transfer between the two sites. 
Then, the frequency of the targeted transfer is 

^A 

" T= 7m (24) 

where K(fc) is the first kind complete elliptic integral. 

ujt is minimum and proportional to the coupling A when xo = 0, i.e., when the 
"nonlinear tuning" between the two sites is optimal. Then the oscillation of the 
difference of action I(t) between the two oscillators becomes a pure cosine. In 
this case, the explicit solution describing the optimal targeted transfer from the 
donor to the acceptor obtained for E = Et = —ujuq/x and Iq — It/2 — —ou/x 
is rather simple 



-UJ ,. 1 . CUo^n 

— cos At exp — 1( — sin 2 At H ) 

X 4A 2 



ty2\t) —i\ sin At exp — 1{— sin 2At H ) 

V X 4A 2 



(25) 
(26) 
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Fig. 3. 3D plot of the rate A of transfer for the action between the two sites of 
the dimer versus E and Iq for Hamiltonian (9) at xo = 0,x = 2, loq = 0, u = — 1. 
Coupling A = 0.01 is weak. 100% transfer occurs at the top of the peak while away 
from the peak, action transfer becomes negligible. 

Fig. 3 shows the 3D plot of the rate of action transfer A = (max i /(t) — 
mm t I(t))/(2I ) between the two sites as a function of the energy E and the 
total action I Q of the solution. 

The targeted energy transfer solution disappears when its frequency reaches 
zero, i.e., k becomes equal to 1. Then, polynomial V(I) gets a degenerate 
zero at I = which means that a time-periodic solution at / = appears. It 
corresponds to a new DB solution which suddenly blocks the targeted transfer. 
At the limit k — 1, the elliptic cosine becomes a sech function [40] and close 
to this limit, the transfer of energy becomes intermittent. 

The time-periodic solutions (DBs and multibreathers) are characterized by 
I = in Eq.(17) and 9 = dH/dl = yields 9 = or tt and 

4A 2 4 2 - ( X oh + Xh + - II) = 0. (27) 
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0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 




Fig. 4. Transfer rate in action versus Iq for x = 2, lo = — 1, A = 0.01 and several 
values of xo = —0.8 (thick dashed), —0.5 (thin dashed-dotted), (thick full line), 
0.25 (thick dashed), 0.5 (thin dashed-dotted), 5. (thin full). 



This is a fourth degree equation which may have 2 or 4 real zeros, functions 
of Iq. The energy of these solutions -Ef>(/o) determines the boundaries of the 
domain in the plane of initial conditions E, I that can be realized. 

We can also analyze partial energy transfer by considering the set of initial 
conditions where all the energy is at site 1. These conditions are characterized 
by 1(0) = I and E = (xo + xKo + ( w + ^oKo- Then, polynomial (19) becomes 

V(I) = \{I - Jo) ((/ - h)(\xo(I + h) + Xh + cu) 2 + 4A 2 (/ + /)) . (28) 

Then, I(t) oscillates between J and the largest zero I m of this polynomial. 
Fig. 4 shows the rate of transfer in action (7 — I m ) I (2^o) versus Iq for a given 
set of parameters XP^i ^ an( i f° r several values of xo- 100% transfer is obtained 
at Iq = It/2 = 1/2 when condition (15) is fulfilled, i.e., \xo\ < 4|Ax|/|a>| = 
0.08. There are discontinuities on these curves where the oscillation regime of 
I(t) becomes cnoidal which are due to the appearance of new time-periodic 
solutions blocking the energy transfer. 
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5 Detecting targeted energy transfer in complex models 

This theory of targeted transfer can be used in practice in complex systems to 
detect targeted transfer because the functions involved in the effective dimer 
Hamiltonian (3) can be explicitly calculated. For that purpose, the initial 
system has to be separated into two coupled subsystems. The first one will 
sustain the donor DB and the second one the acceptor DB. The coupling 
between the two systems is not necessarily weak if the considered DBs are far 
from each other but the effective coupling has to be weak. 

Using the numerical techniques developed earlier for the calculations of DBs 
in arbitrarily complex systems [3-5] , one chooses a DB solution in the isolated 
donor system and calculate its dynamical configuration {pf (u>it) , qf (cuit)} , 
its energy, and its action. Then, varying the frequency, one obtains function 
Hd(I). The same can be done for a DB in the acceptor system. For a complete 
targeted transfer, we require that both DBs be gapless, i.e., Hd(I) and Ha(I) 
can be continued to I = where H D (0) = Ha(0) = 0. 

Since the effective coupling between the two DBs is supposed to be small, It 
can be determined approximately by the equality Ho{It) = Ha(It), which 
requires that the two curves intersect at / 7^ 0. For having targeted transfer 
at reasonably small coupling, function e(I) = H D {I T /2 + I) + Ha(It/2 — I) 
should be small. Once these conditions are fulfilled, we have to check that 
our dimer assumption is valid, i.e., that there is no other DB either on the 
donor system or the acceptor system with almost the energy functions Hd(I) 
or Ha(I) in the interval [0,/r]- 

Next, the coupling energy between the two systems is considered. It is origi- 
nally a function of the global system configuration {pf , qf} and {pf , qf} It 
yields the effective coupling in the dimer model by averaging this energy over 
the fast variable #0 fo r the system configuration {pf ((9q+9)/2), qf \{9q+9)/ '2)} 
and {pf((9 - 0)/2), qf((9 - 9)/2)}. It is thus only a function of I , I and 9 
which has to be small. 

Comparing the precise detuning function and the coupling functions as in 
Eq.(7), we can predict rather accurately the possible existence of targeted 
transfer solutions between these donor and acceptor DBs. In the cases we 
tested up to now, where these conditions are correctly fulfilled, molecular 
dynamics confirmed the existence of targeted transfer where almost 100% of 
the energy oscillates for a very long time between the donor and the acceptor 
DB. As expected also, this energy transfer is blocked as soon as the initial 
amount of energy is not optimal or the coupling too weak [43]. 
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6 Irreversibility in the energy targeted transfer 

DBs are out of equilibrium excitations. In real systems, which are not strictly 
at K, they interact with the phonon modes and other excitations which 
should make their lifetime finite. Thermalization should be in principle a conse- 
quence of intrinsic nonlinearities of the system. However, standard approaches 
usually simulate thermalization and relaxation by a white Langevin noise and 
an extra damping. This is known to be equivalent to a coupling of the initial 
system with a fictitious extra system consisting of an infinite set of harmonic 
oscillators with a uniform distribution of frequencies [44] . This is clearly a bias 
which makes the DBs to systematically radiate energy in this bath and artifi- 
cially washes out their existence as very long lifetime solutions. Actually, the 
long lifetime of DBs is precisely due to the fact that they have no resonance 
with the linear phonon spectrum which implies that either their frequencies 
are above the phonon spectrum or there are gaps in this phonon spectrum. 
It is the nonuniform and specific phonon spectrum of the model itself which 
should be considered for studying DB lifetime as we started to do in some 
recent works [23,28]. In other words, real systems exhibit a selective damping 
depending on the frequencies of their modes which should be explicitly taken 
into account to be physically correct. Thus, the thermal bath assumption with 
a uniform phonon spectrum is not physically founded and, worst, induces the 
denaturation of the model. 

We discussed above the ideal case for targeted energy transfer at K. For 
obtaining an ideal transfer, besides the adiabaticity hypothesis, the donor 
and acceptor DBs must be gapless. This condition requires that there exist no 
linear modes in the donor (or the acceptor) system with frequency in the range 
of variation of the DBs frequency. Indeed, it has been shown in [3-5] that before 
the frequency of a DB collides with a linear mode frequency, it systematically 
bifurcates with another one. In this case, the DB is not strictly continuable to 
zero amplitude. There is a frequency discontinuity corresponding to the region 
where the resonance with the linear mode destroys the DB. Strictly speaking, 
when there is a linear mode in the range of variation of frequency of a DB 
on the donor (or acceptor) system, this DB exhibits a nonvanishing energy 
gap which in principle invalidates our assumptions for describing the coupled 
system as a generalized dimer model (3) . As a consequence, the frequency of 
the ideal targeted transfer solution precisely oscillates in between the edges of 
a gap of the linear phonon spectrum. 

We would like to briefly discuss here what can happen to targeted energy 
transfer if there is just a single resonant linear mode in the phonon gap, first 
in the case of strong coupling, i.e., this linear mode is spatially located close 
to the donor (or the acceptor) DBs. As we just said, the DB of the donor (or 
the acceptor) will exhibit an energy gap, so that the detuning function cannot 
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be defined for the whole interval of variation of I. The energy transfer will 
become incomplete and irreversible. As the frequency of the targeted energy 
transfer solution approaches the linear resonance, a substantial fraction of the 
energy will be suddenly transferred to this linear mode instead of the acceptor 
DB. Consequently, the frequency of the donor DB will change independently of 
the acceptor DB. The loss of resonance between donor and acceptor DB stops 
immediately the energy transfer. If one assumes that the energy captured by 
the linear resonant mode does not return significantly to the donor DB, the 
part of the energy left on the donor system and of the linear mode, could be 
dissipated into phonons through chaotic processes and finally heat. Then, only 
a well-defined fraction of the initial energy on the donor has been transferred 
and the excess of energy has been dissipated into heat. We obtained only an 
incomplete but irreversible targeted transfer. 

The effect should be much less sharp when the linear mode which is resonant 
with the donor DB is weaker, i.e., the resonant linear mode is physically located 
far away from the donor and acceptor DBs, say, in order to fix the ideas, on 
the donor system. Then, there is only a small discontinuity gap in frequency 
to jump for the donor DBs where function Hd(Id) is undefined. We think 
that despite a small energy loss when crossing the resonance, targeted energy 
transfer could be continued beyond this resonance and be almost completed. If 
there are many such weak resonances they could be described as a frequency 
dependent damping depending on the distribution and coupling with reso- 
nant linear modes. Resonances with bands of extended phonons should have 
a similar effect. 

Actually, in realistic systems there are likely some linear modes in the phonon 
gap causing some energy dissipation and making energy transfer irreversible. 
On the opposite, at finite temperature, thermal activation of the linear modes 
could also help targeted energy transfer by stochastic resonance effects. These 
problems will be investigated in the future. 



7 Concluding remarks and new perspectives 

We have demonstrated here that the essential difference between linear and 
nonlinear resonance (allowing the targeted energy transfer) is that nonlinear 
resonance is selective in energy (or action), unlike linear resonance, which is 
not. Targeted energy transfer occurs because the frequencies of the donor os- 
cillator or DB and acceptor oscillator or DB, although they both vary during 
the transfer, remain almost equal. This condition can be achieved at weak cou- 
pling, only for a certain initial energy Et and when a certain detuning function 
is small enough compared to the coupling energy. In the ideal case, an amount 
of energy E T deposited on the donor DB oscillates slowly and periodically 
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back and forth between donor and acceptor. It is totally transferred between 
the donor and the acceptor DB at periodic times. Otherwise, although there 
is some energy oscillation, the energy transfer becomes negligible. 

We do not believe that the obtained targeted energy transfer solutions corre- 
spond to strictly exact solutions of the initial model because of the hypothesis 
of adiabaticity which cannot be perfectly fulfilled. However, situations close 
to the ideal case could exist where the number of energy oscillations between 
donor and acceptor DBs is very large. The high sensitivity to the initial energy 
makes that this energy transfer can be easily stopped by intermediate linear 
resonance which can be viewed as a kind of transistor effect. This effect could 
be used in practical devices (e.g., with optical fibers). 

Within the perspectives opened up by this approach, we would like to suggest 
a series of new physical problems to be investigated in the near future. 

We may consider situations where, instead of one, there are several DBs on 
the donor system with almost identical functions H D {I) and one acceptor 
DB on the acceptor system . Then, if we are in the condition of targeted 
energy transfer for a single donor DB, any donor DBs excited at the energy 
Et will transfer their energy to the acceptor DB. However, the transfer has 
good chances to be irreversible because when the transfer is completed, the 
acceptor energy could return simultaneously to all donor DBs but then the 
resonance between the donor DBs and the acceptor DB will be broken and 
the transfer will be blocked. This situation is interesting for realizing energy 
funneling analogous to what is observed in chlorophyll [45]. 

We may also consider situations such that for an energy E' T < Et-, the accep- 
tor DB (A) is tuned to another acceptor DB (A'). Then, before the transfer 
from the initial donor (D) to the acceptor (A) is completed, this acceptor (A) 
transfer its energy to the next acceptor (A') with some energy loss (note that 
the DB A' acts initially like an intragap linear resonance) but instead of dis- 
sipating energy, it may collect most of the A acceptor energy. By this way, 
recurrent transfer could be built along a certain well determined sequence of 
acceptor DBs. We could then realize irreversible cascades of energy transfer 
(with some energy losses however). It is also important to realize now that 
DBs in molecules may have piezoactive effects, i.e., they might generate de- 
formation of the molecules as emphasized in [20] . 

This effect could be important for example for understanding the energy con- 
version of ATP into mechanical energy (molecular motors). The energy depo- 
sition from ATP could occur first as a DB located at a certain receptive site of 
the moving biomolecule (e.g., kinesin adsorbed on a microtubule). Then, this 
DB could cascade by targeted energy transfer through a series of DBs located 
at specific sites of the biomolecule generating a well determined sequence of de- 
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formations which lead to a well determined motion of the biomolecule (step) 
along the oriented microtubule after the time t m necessary for this cas- 
cade process. Beyond this time, the DB energy is finally relaxed into heat. 
Thus, excitations by ATP randomly repeated at time intervals longer than 
this characteristic time t m should produce systematically a directive walk of 
this biomolecule along the microtubule (this motion should become randomly 
intermittent in time if the ATP concentration is low). In the future, we plan to 
use and to develop our tools for calculating DBs in nonlinear systems for con- 
trolling precisely their targeted energy transfer, funneling, cascades, piezoac- 
tivity of DBs and thus for designing artificial models with desired properties 
mimicking those of biosystems. 

Let us also note that instead of being atomic vibrations, these DBs could 
be electronic vibrations (excitons), polarons, or atomic vibrations combined 
with an electronic charge (polarobreathers) [17]. Thus, there are many possi- 
ble variations of the targeted transfer principle depending on the considered 
physical problem (transfer of energy, charge, spin, etc) and its scale of energy. 
This theory of targeted transfer could be useful, for example, for revisiting the 
standard theories of fluorescence ([46]), chemical reactions, catalysis,etc. 
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